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A conserved current for generic quadratic curvature gravitational models is denned and it is shown 
that at the linearized level it corresponds to the Deser-Tekin charges. An explicit expression for 
the charge for new massive gravity in three dimensions is given. Some implications of the linearized 
equations are discussed. 
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I. INTRODUCTION 



Defining a local energy density for a given gravitating system is still a matter of discussion, however, it is possible 
to define a total energy in the general theory of relativity [1]. Closely following the definition of conserved charges 
in the general theory of relativity [2], a new definition of energy for the quadratic curvature (QC) models has been 
proposed in [3] which will be called Deser-Tekin (DT) charges (see, also the previous definitions [4]). In particular, 
the conserved charge corresponding to the DT energy for asymptotically flat or constant curvature geometries have 
desirable properties, for example, that it circumvents the zero-energy theorem for conformal gravity in four dimensions 
[5]. 

Sparling forms [6] are also used to define energy in the general relativity which are well-suited to a treatment in terms 
of differential forms relative to an orthonormal coframc [7] . In a direction parallel to the present work, (generalized) 
Sparling forms for dimensionally-continued Eulcr forms in arbitrary dimensions were introduced in [8]. The Sparling 
forms can also be related to the so-called Arnowit-Deser-Misner (ADM) energy [1] relative to a coordinate basis and 
other energy- momentum pseudo-tensors such as Landau-Lifschitz energy pseudo-tensor as well [9]. They also play 
an important role in the hamiltonian formulation of General theory of relativity [10] and even in the analysis of 
gravitational radiation in post- Newtonian approximation for theories with torsion [11]. In the present work, by using 
the language of the exterior algebra and orthonormal coframe formulation for the QC gravity models, a conserved 
quantity will be defined and its connection with the DT energy will be discussed. 

The rest of the paper is organized as follows. In the following section, metric QC field equations are described 
in a form that is useful in defining an on-shell conserved current and a total energy. In the subsequent section, the 
linearized QC equations are studied to further investigate the relation of the energy definition presented with the DT 
charges. In the last section, an explicit expression for the conserved current in the form of an exact differential form 
for New massive gravity in three dimensions is presented. 



II. QC FIELD EQUATIONS AND DEFINITION OF A NEW CONSERVED CURRENT 

The notation and the conventions for the geometrical quantities closely follow those of [12, 13], see also [9]. For 
the technical details related, in particular, to the QC field equations relative to an orthonormal coframc, the reader 
is referred to [12]. 

A convenient starting point is to recall that the identity *G a = —^Q^ A *Q a ^ u for the Einstein form can be used 
to split it into the sum of two terms as 

*G a =d*F a + *t a (1) 
where F a is known as Sparling 2-form [6] with the expression 

*F a = A *0°*"' (2) 
whereas the pseudo-tensorial energy-momentum 3-form *t a is explicitly given by 

*t a = §(<«v A u a x A *6> xvX - w„ A A u\ A *9^ a ). (3) 
The standard ADM energy expression for an asymptotically flat spacetime follows from the integral of *F°, namely 

where £ is a spacelike hypersurface (see, e.g., [9], p. 118). For the well-known example of Schwarzchild metric, by in- 
troducing the isotropic coordinates, the formulae (4), singling out the leading order terms in the metric asymptotically, 
yields the mass parameter of the metric [9] . 

As advertised above, it is possible to define conserved current, or equivalently an exact form, for generic QC gravity 
by splitting the field equations in a way similar to (1), cf. Eqn. (13) below. For the sake of clarity for the arguments 
and the definitions, mainly pure QC gravity models that follow from the general lagrangian density of the form 

C qc = aR 2 *l + bR a A *R a (5) 

will be considered (a, b are dimensionless coupling constants), confining the discussion to four and subsequently to 
three dimensions. For the use of first order formalism with the constraints Lo a /3 + up a — (metric compatibility) and 
a = (zero-torsion), one starts with the extended lagrangian 

C e = C qc + A" A 6 Q (6) 
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where A Q is a vector- valued 2-form imposing the dynamical constraint G Q = 0. The metric compatibility constraint 
can simply be incorporated into the variational procedure by 5uj a /3 +6LU/3 a — 0. Subsequently, by taking the variational 
derivative of C e [9 a , , X a ] , one obtains metric equations from the coframe equations expressed in terms of the pseudo- 
Ricmannian quantities only. For the details of the derivation of the QC field equations relative to an orthonormal 
coframe refer to [12, 13]. It follows that the vacuum equations for the QC lagrangian density (5) then can be written 
in the form 



*E a = DX a + *T a = (7) 

where the vector valued 1-form E a — E a ^9^ is defined by (7) and is a symmetric vector-valued 1-form (by invariance 
of the action under local coframe rotation). As in the Einstein-Hilbert case, it is covariantly constant, D * E a = 
as a result of diffeomorphism invariance of the QC lagrangian. Moreover, the pure QC terms are contained in *T a 
term where T a = is another symmetric vector-valued 1-form. In a generic QC model, fourth order lagrange 

multiplier term, DX a , and second order terms, *T Q , are both expressible in terms of auxiliary antisymmetric-tensor- 
valued 2-form X al3 . For the QC lagrangian (5) it is explicitly given by 

The connection equations can be solved uniquely for lagrange multiplier and consequently, the lagrange multiplier 
2-forms A' 3 can be found in terms of the auxiliary 2-form as 

A* 3 = 2i a U aP + \6 f3 A i^i v W v 

= 2 * D[bR + (2a + ±b)R9 fj ] (9) 

with H a/3 = D * X a P whereas the second term in (7) has the explicit form 

*r = A i a * X„ v - i a W A *X IM/ ). (10) 

Note that the trace of *T a is T a a * 1 = (4 — n)C qc in n dimensions. The form of the field equation (7), which has 
been introduced in [12], directly allows one to define a conserved quantity by inspection. In order to do so, first note 
that the fourth order term can simply be written out as 

DX a = dX a + uj a p A X 13 (11) 

which in fact separates the terms containing fourth order partial derivatives of the metric components from the 
lower orders relative to a coordinate basis. This then motivates the definition of a convenient vector-valued 1-form 

r a = T a 6P by 

*r a = *T a + Lu a p A X fj . (12) 

Thus, by making use of the definitions (11) and (12) in the vacuum QC gravity equations, (7) then leads to the 
compact equations as 



dX a + *r a = 0. (13) 

Consequently, the straightforward splitting renders *T a an exact form and, by the operator identity d 2 = 0, it is also a 
closed form d*r a = 0. Note however that *r Q is , as in the case of the Sparling form, pseudo-tensorial since the right 
hand side of (12) depends explicitly on the connection 1-form. On the other hand, for QC gravity, A" will always be 
linear in curvature components. Therefore, the energy-momentum 3-form (12) qualifies to define a conserved current 
for generic QC models (5). Consequently, it is possible to define the four momentum P a as 

P Q (E) = [ X a = — f *t q (14) 

where S is a space-like hypersurface with boundary 9S. By definition, P Q (S) vanishes identically for dT. ~ 0. The 
definition of *r" is independent of matter content in the sense that any matter energy-momentum can simply be 
included in the right hand side of (12) provided that it is covariantly constant. By inspecting the QC field equations 
written out relative to a coordinate basis (see, e.g., Eqn. (3) in Ref. [3])), one concludes that it is not straightforward 
to obtain the splitting (13) relative to a coordinate basis. 
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III. LINEARIZED QC EQUATIONS 

A remarkable property of the conserved current A Q appears to be related to the linearized form of the vacuum 
equations (7). In order to establish connection with explicit expression defining energy presented in [3], let us consider 
the linearization of the field equations (7) by introducing the perturbation h aj 3 to the metric g — rj a [j9 a ® 0^ with 
Va/3 = diag( — h ++)• In terms of a local orthonormal basis of coframe 1-forms, they can be defined as 

e a Kdx a + h a dx fl . (15) 

The indices of the symmetric h a p are raised and lowered by the flat Minkowski metric r\ in the linear approximation. 
The relevant tensorial quantities which are linear in h a p below are indicated by a label L. The harmonic gauge, where 
the coordinates at use satisfy d * dx a — 0, leads to the equation d a h a fj — ^dph"^ = and the linearized curvature 
forms in this gauge are 

R L a = -Eih a pdx p , R L = -nh a a (16) 

where □ = rj al3 d a dj3 for convenience [14]. From these expressions, it follows that the linearized Einstein 1-form is 
G a = R a - \R L dx a or equivalently *G% = d* F£. The contracted Bianchi identity D * G a = holds at the 

linearized level in the form d*G1 = dpG^ * 1 — 0. With these preliminary considerations, the vacuum QC equations 
(7) can now easily be linearized around Minkowski background. (It is possible to linearize (7) around a constant 
curvature background which is also known to solve (7) exactly, see, e.g. [12]. Constant curvature backgrounds will 
not be considered in the discussion here) Since *T a term involves only QC contributions, the linearization in this case 
involves only the constraint term DX a , and writing out this term explicitly one finds 

*E a = 2D * D[bR a + (2a + \b)R6 a ] + *T a = 0. (17) 

Therefore, after evaluating the exterior covariant derivatives in the constraint term and writing all the terms out 
explicitly, *E a = then takes the form 

2d * d[bR a + {2a + \b)R6 a ] + 2d * {uj a A [bR 13 + {2a + ±b)R6 p }} 

+ 2w a p A *d[bR f3 + (2a + ±b)R6 fj ] + 2u a fj A *{w /S „ A [bR" + (2a + ±b)R6>*}} 

+ n» v Ai a * [9^ A (bR u + aR6 u ) - 6 V A (bR^ + aR6^)] - i a C qc [a, b] = 0. (18) 

DT charges were defined by expanding this equation around a constant curvature or flat background. In particular, 
the DT energy is obtained by linearizing the (00) component of the equation. As will explicitly be shown below, the 
definition of DT charge follows from the linearization of the first term in (18). Only the first term in (18) contains 
the terms linear in the partial derivatives of the metric components and it is remarkable that this particular term is 
an exact differential form in the full QC field equations. An important observation is that the linearization of the 
first term is obtained from DX a simply by replacing the covariant exterior derivatives with exterior derivatives and 
by adopting the linearized curvatures and the Hodge dual. Futhermore, by approximating the orthonormal coframe 
1-forms as a ~ dx a , the QC equations linearized in metric components then can be written as 

*E<£ = d\ a L = 2d * d[bR1 p dx fj + (2a + \b)R L dx a ] . (19) 

In (19) and in all the linearized equations involving Hodge dual below, it is necessary that the Hodge dual is also 
to be approximated by that of the flat Minkowski space, for example, 9 a A w dx a A *dx 13 = r] al3 * 1, etc. Now 
assuming that (19) is an expression in the flat background, then (19) explicitly yields 

*E2 - 2d v d^bR a Lp + (2a + \b)5^R L ]dx v A *(dx» A dx?). (20) 

By evaluating the inner products of the basis 1-forms indicated on the right hand side of (20) with the help of the 
identity [9] 

dx v A *(dx" A dx ) = -jy"" * dx p + rf !i * dx" (21) 

one eventually ends up with 



*EZ = 2P a »[bG L ^ + (2a + b)^R L ] * dx fi 



(22) 
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where, following the notation in [3], the linear differential operator P a p with respect to Minkowski background defined 

by 

Pa/3 = d a d/3 - f] a0 D (23) 

is introduced for convenience. The crucial expression (23) for the projection operator P a p results essentially from the 
contractions in the identity (21) which in turn follows from the general expression of A Q in (9). In the above notation, 
P a p can equally be defined by the action of (flat space) operator d * d on a vector- valued 1-form a a = a a j3 dxP by 
d * d<r a = *a' a with a' a = P^pa al3 dx^ . It follows from the definition P a p and linearized Bianchi identity d a G1g = 
that PapG^p = —OG^p and therefore dropping the Hodge duals on both sides to write (22) component form, one 
obtains 

E L afj = -2\bUG% + {2a + b){ Vafj U - d a dp)R L }. (24) 

The expression in the square brackets is precisely the left hand side of Eqn. (5) in [3]. As it was shown in [3], one of the 
crucial properties of the (00) component of the expression is that it involves only the second order time derivatives of 
metric as a result of the fact that Poo = — V 2 and consequently it corresponds to a Poisson type potential. Therefore 
the DT energy is related to the asymptotical behavior curvature rather than the asymptotical behavior of metric. 
Under certain assumptions about spatial asymptotics, (14) is therefore expected to yield results identical to DT 
charge by singling out leading terms of the curvature. Finding a DT charge for a given (asymptotically flat) solution 
of the QC equations requires the curvature to vanish as ~ r _1 at spatial infinity and consequently, such solutions 
exclude the metrics having the asymptotic behavior ~ r _1 [3]. Thus, one concludes that, without resorting to linear 
approximation, it is possible to define a conserved current for QC gravity and the DT energy definition corresponds 
to the linearized Lagrange multiplier term d\° L under certain assumptions about asymptotic behavior. Finally, note 
that the linearization of the current defined in (12) then corresponds to side of Eqn. (5) in [3] up to a Hodge 
dual. 

Now a brief discussion about the consistency of linearized equations and some properties of their solution is in 
order. Any solution of vacuum Einstein field equations is also a solution of the vacuum QC field equations (7), that 
is, DX a = and *T a — separately because of X a/3 — identically for a metric with R a — in the above notation. 
However, the exact solutions of fourth order gravity has more integration constants [15] and even the linearized 
equations (and the solutions to them) have peculiar properties [16]. 

The property that the field equations are covariantly constant, D * E a = 0, continues to hold in the linear approx- 
imation. Explicitly, the linearized Bianchi identity reads 

d a E a Lp = 2d a P a » [bG L ^ + (2a + fo)^P L ] (25) 

where d a E1 B = identically follows from the operator identity d a P a ^ = by definition (23). In the presence of 
matter field sources, the field equations *E a = — 2 * T™ coupled with matter energy-momentum 1-form T™ = T™ g 9^ 
requires that at the linearized level, d a T™g = 0. The consistency of the linearized equations (22) require a further 
differential constraint on matter energy-momentum T™ 3 (which is not present in the full non-linear equations). By 
calculating the trace of (22) and making use of the definition (23) one finds that linearization of *E a = — 2 * T™ 
implies 

DT "™ + Wa^b) PapTm = hU2G ^ (26) 

where T m stands for the trace of the matter energy-momentum tensor. For the subcase b = 0, the constraint (26) 
was derived long ago [15] to show that linearized equations for R 2 gravity do not support solutions with the energy- 
momentum tensor of the form T™p — <5°<5°p(x) for static, extended (and bounded) matter distribution p(x). Explicitly, 
for the subcase 6 = 0, (00) component of (26) implies p(x) = [15]. For the general case 6^0 however, it yields 
p(x) ~ AGqo where P o = rfididj = A and the Latin indices run over flat spatial coordinates. This is in accordance 
with the fact that the expression inside the square bracket in (22) can be regarded as a Poisson potential for the linear 
approximation of the theory [3]. 



IV. A CONSERVED CURRENT FOR NEW MASSIVE GRAVITY 

The QC field equations (7) have the same form in any dimensions n > 3. For a — n/4(n — 1) and b = 1 in (5), the 
expression for the lagrange multiplier form in (9) becomes proportional to Cotton 2-form. In particular, for n = 3, 
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the lagrangian density for New massive gravity (NMG)[17] reads 

Cnmg = -R*l + \Ck (27) 

up to an overall coupling constant which has been dropped for convenience. The parameter m has necessarily the 
dimension of mass and the QC part is Ck — C qc [a — — |, b = 1]. At the linearized level (27) is equivalent to Pauli-Fierz 
lagrangian for a massive spin-2 particle. 

In accordance with the above description, it is convenient to write the field equations for the QC part of (27) in 
the form given in (7) as *E^- = 0. The vacuum equations that follow from (27) then take the form 

Each of the terms in (28) can be split to define an exact pseudo-tensorial 2-forms, for the QC part and *t a that 
is given in (3) for Einstcin-Hilbert part. Consequently, (28) can be rewritten in a Maxwell-like form as 

d * (F a + m- 2 C a ) + *r£ MG - (29) 

where *r^ MG explicitly reads 

*r£ MG - *t a + ^ Kj A *C" + \ * T£) (30) 

by making use of the fact that A Q = 2m~ 2 * C a with C a standing for the Cotton 2-form. On the right hand side of 
(30), *T£ = *T a [fl^,X^} with X£ v = X» v [a = -§,& = 1]. Consistently, at the linearized level, Eqn. (29) reads 

d* (FjZ +m- 2 C2) = (31) 

(also assuming that at the linearized level * stands for Minkowski Hodge dual) and (31) yields the field equation for 
a propagating massive spin-2 particle [17]. 

It is worth to emphasize that in (29) the *F° component of the Sparling form (the asymptotically dominant term) 
is related to ADM energy whereas A = mT 2 * C° component is related to DT energy definition for generic higher 
curvature gravity. The new energy definition in Section II allows one to combine these energy definitions for NMG in 
a natural way. 

As mentioned in the introduction, one of the motivations for the construction of the DT charges is to solve the zero 
energy problem of Weyl gravity with flat background. Weyl gravity follows from the lagrangian C qc [a = —\,b = 1] 
in four dimensions and in this case, \ a = 2 * C a as in NMG [12]. Consequently, in four dimensions P a in (14) does 
not vanish identically in accordance with the DT result. On the other hand, for n > 4 dimensions, A Q = identically 
for Gauss-Bonnet lagrangian density n a p A VL^ V A *9 af>tJ - v and the corresponding field equations are second order in 
metric components. Consequently, in this case P° in (14) yields zero identically for the energy as the DT definition 
of energy does as well. 

The author would like to thank Ozgur Delice for helpful comments and critical reading of the manuscript. 
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